Search for: All records

Abstract Insect wings are heterogeneous structures, with flexural rigidity varying one to two orders of magnitude over the wing surface. This heterogeneity influences the deformation the flapping wing experiences during flight. However, it is not well understood how this flexural rigidity gradient affects wing performance. Here, we develop a simplified 2D model of a flapping wing as a pitching, plunging airfoil using the assumed mode method and unsteady vortex lattice method to model the structural and fluid dynamics, respectively. We conduct parameter studies to explore how variable flexural rigidity affects mean lift production, power consumption and the forces required to flap the wing. We find that there is an optimal flexural rigidity distribution that maximizes lift production; this distribution generally corresponds to a 3:1 ratio between the wing’s flapping and natural frequencies, though the ratio is sensitive to flapping kinematics. For hovering flight, the optimized flexible wing produces 20% more lift and requires 15% less power compared to a rigid wing but needs 20% higher forces to flap. Even when flapping kinematics deviate from those observed during hover, the flexible wing outperforms the rigid wing in terms of aerodynamic force generation and power across a wide range of flexural rigidity gradients. Peak force requirements and power consumption are inversely proportional with respect to flexural rigidity gradient, which may present a trade-off between insect muscle size and energy storage requirements. The model developed in this work can be used to efficiently investigate other spatially variant morphological or material wing features moving forward. 

Flapping insect wings experience appreciable deformation due to aerodynamic and inertial forces. This deformation is believed to benefit the insect’s aerodynamic force production as well as energetic efficiency. However, the fluid-structure interaction (FSI) models used to estimate wing deformations are often computationally demanding and are therefore challenged by parametric studies. Here, we develop a simple FSI model of a flapping wing idealized as a two-dimensional pitching-plunging airfoil. Using the Lagrangian formulation, we derive the reduced-order structural framework governing wing’s elastic deformation. We consider two fluid models: quasi-steady Deformable Blade Element Theory (DBET) and Unsteady Vortex Lattice Method (UVLM). DBET is computationally economical but does not provide insight into the flow structure surrounding the wing, whereas UVLM approximates flows but requires more time to solve. For simple flapping kinematics, DBET and UVLM produce similar estimates of the aerodynamic force normal to the surface of a rigid wing. More importantly, when the wing is permitted to deform, DBET and UVLM agree well in predicting wingtip deflection and aerodynamic normal force. The most notable difference between the model predictions is a roughly 20° phase difference in normal force. DBET estimates wing deformation and force production approximately 15 times faster than UVLM for the parameters considered, and both models solve in under a minute when considering 15 flapping periods. Moving forward, we will benchmark both low-order models with respect to high fidelity computational fluid dynamics coupled to finite element analysis, and assess the agreement between DBET and UVLM over a broader range of flapping kinematics.